Examination of the Potential of Terrestrial Laser Scanning and Structure-from-Motion Photogrammetry for Rapid Nondestructive Field Measurement of Grass Biomass

Above ground biomass (AGB) is a parameter commonly used for assessment of grassland systems. Destructive AGB measurements, although accurate, are time consuming and are not easily undertaken on a repeat basis or over large areas. Structure-from-Motion (SfM) photogrammetry and Terrestrial Laser Scanning (TLS) are two technologies that have the potential to yield precise 3D structural measurements of vegetation quite rapidly. Recent advances have led to the successful application of TLS and SfM in woody biomass estimation, but application in natural grassland systems remains largely untested. The potential of these techniques for AGB estimation is examined considering 11 grass plots with a range of biomass in South Dakota, USA. Volume metrics extracted from the TLS and SfM 3D point clouds, and also conventional disc pasture meter settling heights, were compared to destructively harvested AGB total (grass and litter) and AGB grass plot measurements. Although the disc pasture meter was the most rapid method, it was less effective in AGB estimation (AGBgrass r2 = 0.42, AGBtotal r2 = 0.32) than the TLS (AGBgrass r2 = 0.46, AGBtotal r2 = 0.57) or SfM (AGBgrass r2 = 0.54, AGBtotal r2 = 0.72) which both demonstrated their utility for rapid AGB estimation of grass systems

Comparisons of Temporal and Spatial Trends in the Spatially Complete Global Spectral Surface Albedos Products

Five years of spatially complete snow-free land surface albedo data have been prepared using high quality white-sky and black-sky land surface albedo observations (MOD43B3) from the MODIS instrument aboard NASA's Terra satellite platform. The data were generated using an updated ecosystem-dependent temporal interpolation technique. In this paper, we describe the refinements in the technique and the creation of a spatially complete snow-free five-year aggregate climatology product. This paper also describes an error analysis of the interpolation technique. The filled albedo products are examined through comparisons of temporal and spatial trends for pixels that have been filled versus pixels that have been retained in the original MOD43B3 values. The variability in the trends showcase how the filling technique maintains the pixel-level spatial, spectral, and temporal integrity of the MOD43B3 data. These comparisons are made for both a single-year of filled data, year 2002, and for the five-year aggregate climatology product

Quantifying the effects of freeze-thaw transitions and snowpack melt on land surface albedo and energy exchange over Alaska and Western Canada

Variations in land surface albedo and snow-cover strongly impact the global biosphere, particularly through the snow-albedo feedback on climate. The seasonal freeze-thaw (FT) transition is coupled with snowpack melt dynamics and strongly impacts surface water mobility and the energy budget in the northern (≥45°N) arctic and boreal region (ABR). However, understanding of the regional variation in snowmelt and its effect on the surface energy budget are limited due to sparse in situ measurements of these processes and environmental constraints on effective monitoring within the ABR. In this study, we combined synergistic observations from overlapping satellite optical-infrared and microwave sensor records to quantify the regional patterns and seasonal progression in wet snow conditions during the spring snowmelt and autumn snow accumulation periods across Alaska and western Canada. The integrated satellite record included daily landscape FT status from AMSR microwave brightness temperature retrievals; and snow-cover extent, black sky albedo and net shortwave solar radiation (R snet) derived from MODIS and AVHRR observations. The integrated satellite records were analyzed with in situ surface air temperature and humidity observations from regional weather stations over a two-year study period (2015–2016) overlapping with the NASA ABoVE (Arctic Boreal Vulnerability Experiment). Our results show a large (79%) mean decline in land surface albedo between dry snow and snow-free conditions during the spring (March–June) and autumn (August–November) transition periods. Onset of diurnal thawing and refreezing of the surface snow layer and associated wet snow conditions in spring contributed to an approximate 25% decrease in snow cover albedo that extended over a seven to 21 week snowpack depletion period. The lower wet snow albedo enhances R snet by approximately 74% (9–10 MJ m−2 d−1) relative to dry snow conditions, reinforcing snowmelt and surface warming, and contributing to growing season onset and activation of biological and hydrological processes in the ABR. These results contribute to better understanding of snow albedo feedbacks to Arctic amplification, and the representation of these processes in global Earth system models

An Approach for the Long-Term 30-m Land Surface Snow-Free Albedo Retrieval from Historic Landsat Surface Reflectance and MODIS-based A Priori Anisotropy Knowledge

Land surface albedo has been recognized by the Global Terrestrial Observing System (GTOS) as an essential climate variable crucial for accurate modeling and monitoring of the Earth's radiative budget. While global climate studies can leverage albedo datasets from MODIS, VIIRS, and other coarse-resolution sensors, many applications in heterogeneous environments can benefit from higher-resolution albedo products derived from Landsat. We previously developed a "MODIS-concurrent" approach for the 30-meter albedo estimation which relied on combining post-2000 Landsat data with MODIS Bidirectional Reflectance Distribution Function (BRDF) information. Here we present a "pre-MODIS era" approach to extend 30-m surface albedo generation in time back to the 1980s, through an a priori anisotropy Look-Up Table (LUT) built up from the high quality MCD43A BRDF estimates over representative homogenous regions. Each entry in the LUT reflects a unique combination of land cover, seasonality, terrain information, disturbance age and type, and Landsat optical spectral bands. An initial conceptual LUT was created for the Pacific Northwest (PNW) of the United States and provides BRDF shapes estimated from MODIS observations for undisturbed and disturbed surface types (including recovery trajectories of burned areas and non-fire disturbances). By accepting the assumption of a generally invariant BRDF shape for similar land surface structures as a priori information, spectral white-sky and black-sky albedos are derived through albedo-to-nadir reflectance ratios as a bridge between the Landsat and MODIS scale. A further narrow-to-broadband conversion based on radiative transfer simulations is adopted to produce broadband albedos at visible, near infrared, and shortwave regimes.We evaluate the accuracy of resultant Landsat albedo using available field measurements at forested AmeriFlux stations in the PNW region, and examine the consistency of the surface albedo generated by this approach respectively with that from the "concurrent" approach and the coincident MODIS operational surface albedo products. Using the tower measurements as reference, the derived Landsat 30-m snow-free shortwave broadband albedo yields an absolute accuracy of 0.02 with a root mean square error less than 0.016 and a bias of no more than 0.007. A further cross-comparison over individual scenes shows that the retrieved white sky shortwave albedo from the "pre-MODIS era" LUT approach is highly consistent (R(exp 2) = 0.988, the scene-averaged low RMSE = 0.009 and bias = 0.005) with that generated by the earlier "concurrent" approach. The Landsat albedo also exhibits more detailed landscape texture and a wider dynamic range of albedo values than the coincident 500-m MODIS operational products (MCD43A3), especially in the heterogeneous regions. Collectively, the "pre-MODIS" LUT and "concurrent" approaches provide a practical way to retrieve long-term Landsat albedo from the historic Landsat archives as far back as the 1980s, as well as the current Landsat-8 mission, and thus support investigations into the evolution of the albedo of terrestrial biomes at fine resolution

Evaluation of Surface and Near-Surface Melt Characteristics on the Greenland Ice Sheet using MODIS and QuikSCAT Data

The Greenland Ice Sheet has been the focus of much attention recently because of increasing melt in response to regional climate warming. To improve our ability to measure surface melt, we use remote-sensing data products to study surface and near-surface melt characteristics of the Greenland Ice Sheet for the 2007 melt season when record melt extent and runoff occurred. Moderate Resolution Imaging Spectroradiometer (MODIS) daily land-surface temperature (LST), MODIS daily snow albedo, and a special diurnal melt product derived from QuikSCAT (QS) scatterometer data, are all effective in measuring the evolution of melt on the ice sheet. These daily products, produced from different parts of the electromagnetic spectrum, are sensitive to different geophysical features, though QS- and MODIS-derived melt generally show excellent correspondence when surface melt is present on the ice sheet. Values derived from the daily MODIS snow albedo product drop in response to melt, and change with apparent grain-size changes. For the 2007 melt season, the QS and MODIS LST products detect 862,769 square kilometers and 766,184 square kilometers of melt, respectively. The QS product detects about 11% greater melt extent than is detected by the MODIS LST product probably because QS is more sensitive to surface melt, and can detect subsurface melt. The consistency of the response of the different products demonstrates unequivocally that physically-meaningful melt/freeze boundaries can be detected. We have demonstrated that these products, used together, can improve the precision in mapping surface and near-surface melt extent on the Greenland Ice Sheet

Variability in Surface BRDF at Different Spatial Scales (30 m-500 m) Over a Mixed Agricultural Landscape as Retrieved from Airborne and Satellite Spectral Measurements

Over the past decade, the role of multiangle remote sensing has been central to the development of algorithms for the retrieval of global land surface properties including models of the bidirectional reflectance distribution function (BRDF), albedo, land cover/dynamics, burned area extent, as well as other key surface biophysical quantities represented by the anisotropic reflectance characteristics of vegetation. In this study, a new retrieval strategy for fine-to-moderate resolution multiangle observations was developed, based on the operational sequence used to retrieve the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5 reflectance and BRDF/albedo products. The algorithm makes use of a semiempirical kernel-driven bidirectional reflectance model to provide estimates of intrinsic albedo (i.e., directional-hemispherical reflectance and bihemispherical reflectance), model parameters describing the BRDF, and extensive quality assurance information. The new retrieval strategy was applied to NASA's Cloud Absorption Radiometer (CAR) data acquired during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC) over the well-instrumented Atmospheric Radiation Measurement Program (ARM) Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site in Oklahoma, USA. For the case analyzed, we obtained approx.1.6 million individual surface bidirectional reflectance factor (BRF) retrievals, from nadir to 75 off-nadir, and at spatial resolutions ranging from 3 m - 500 m. This unique dataset was used to examine the interaction of the spatial and angular characteristics of a mixed agricultural landscape; and provided the basis for detailed assessments of: (1) the use of a priori knowledge in kernel-driven BRDF model inversions; (2) the interaction between surface reflectance anisotropy and instrument spatial resolution; and (3) the uncertain ties that arise when sub-pixel differences in the BRDF are aggregated to a moderate resolution satellite pixel. Results offer empirical evidence concerning the influence of scale and spatial heterogeneity in kernel-driven BRDF models; providing potential new insights into the behavior and characteristics of different surface radiative properties related to land/use cover change and vegetation structure

Examination of the Potential of Terrestrial Laser Scanning and Structure-from-Motion Photogrammetry for Rapid Nondestructive Field Measurement of Grass Biomass

Above ground biomass (AGB) is a parameter commonly used for assessment of grassland systems. Destructive AGB measurements, although accurate, are time consuming and are not easily undertaken on a repeat basis or over large areas. Structure-from-Motion (SfM) photogrammetry and Terrestrial Laser Scanning (TLS) are two technologies that have the potential to yield precise 3D structural measurements of vegetation quite rapidly. Recent advances have led to the successful application of TLS and SfM in woody biomass estimation, but application in natural grassland systems remains largely untested. The potential of these techniques for AGB estimation is examined considering 11 grass plots with a range of biomass in South Dakota, USA. Volume metrics extracted from the TLS and SfM 3D point clouds, and also conventional disc pasture meter settling heights, were compared to destructively harvested AGB total (grass and litter) and AGB grass plot measurements. Although the disc pasture meter was the most rapid method, it was less effective in AGB estimation (AGBgrass r 2 = 0.42, AGBtotal r 2 = 0.32) than the TLS (AGBgrass r 2 = 0.46, AGBtotal r 2 = 0.57) or SfM (AGBgrass r 2 = 0.54, AGBtotal r 2 = 0.72) which both demonstrated their utility for rapid AGB estimation of grass systems

Use of In Situ and Airborne Multiangle Data to Assess MODIS- and Landsat-based Estimates of Surface Albedo

The quantification of uncertainty of global surface albedo data and products is a critical part of producing complete, physically consistent, and decadal land property data records for studying ecosystem change. A current challenge in validating satellite retrievals of surface albedo is the ability to overcome the spatial scaling errors that can contribute on the order of 20% disagreement between satellite and field-measured values. Here, we present the results from an uncertain ty analysis of MODerate Resolution Imaging Spectroradiometer (MODIS) and Landsat albedo retrievals, based on collocated comparisons with tower and airborne multi-angular measurements collected at the Atmospheric Radiation Measurement Program s (ARM) Cloud and Radiation Testbed (CART) site during the 2007 Cloud and Land Surface Interaction Campaign (CLAS33 IC 07). Using standard error propagation techniques, airborne measurements obtained by NASA s Cloud Absorption Radiometer (CAR) were used to quantify the uncertainties associated with MODIS and Landsat albedos across a broad range of mixed vegetation and structural types. Initial focus was on evaluating inter-sensor consistency through assessments of temporal stability, as well as examining the overall performance of satellite-derived albedos obtained at all diurnal solar zenith angles. In general, the accuracy of the MODIS and Landsat albedos remained under a 10% margin of error in the SW(0.3 - 5.0 m) domain. However, results reveal a high degree of variability in the RMSE (root mean square error) and bias of albedos in both the visible (0.3 - 0.7 m) and near-infrared (0.3 - 5.0 m) broadband channels; where, in some cases, retrieval uncertainties were found to be in excess of 20%. For the period of CLASIC 07, the primary factors that contributed to uncertainties in the satellite-derived albedo values include: (1) the assumption of temporal stability in the retrieval of 500 m MODIS BRDF values over extended periods of cloud-contaminated observations; and (2) the assumption of spatial 45 and structural uniformity at the Landsat (30 m) pixel scale

A Cross Comparison of Spatiotemporally Enhanced Springtime Phenological Measurements From Satellites and Ground in a Northern U.S. Mixed Forest

Cross comparison of satellite-derived land surface phenology (LSP) and ground measurements is useful to ensure the relevance of detected seasonal vegetation change to the underlying biophysical processes. While standard 16-day and 250-m Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation index (VI)-based springtime LSP has been evaluated in previous studies, it remains unclear whether LSP with enhanced temporal and spatial resolutions can capture additional details of ground phenology. In this paper, we compared LSP derived from 500-m daily MODIS and 30-m MODIS-Landsat fused VI data with landscape phenology (LP) in a northern U.S. mixed forest. LP was previously developed from intensively observed deciduous and coniferous tree phenology using an upscaling approach. Results showed that daily MODIS-based LSP consistently estimated greenup onset dates at the study area (625 m × 625 m) level with 4.48 days of mean absolute error (MAE), slightly better than that of using 16-day standard VI (4.63 days MAE). For the observed study areas, the time series with increased number of observations confirmed that post-bud burst deciduous tree phenology contributes the most to vegetation reflectance change. Moreover, fused VI time series demonstrated closer correspondences with LP at the community level (0.1-20 ha) than using MODIS alone at the study area level (390 ha). The fused LSP captured greenup onset dates for respective forest communities of varied sizes and compositions with four days of the overall MAE. This study supports further use of spatiotemporally enhanced LSP for more precise phenological monitoring

Forest structure and aboveground biomass in the southwestern United States from MODIS and MISR

Red band bidirectional reflectance factor data from the NASA MODerate resolution Imaging Spectroradiometer (MODIS) acquired over the southwestern United States were interpreted through a simple geometric–optical (GO) canopy reflectance model to provide maps of fractional crown cover (dimensionless), mean canopy height (m), and aboveground woody biomass (Mg ha−1) on a 250 m grid. Model adjustment was performed after dynamic injection of a background contribution predicted via the kernel weights of a bidirectional reflectance distribution function (BRDF) model. Accuracy was assessed with respect to similar maps obtained with data from the NASA Multiangle Imaging Spectroradiometer (MISR) and to contemporaneous US Forest Service (USFS) maps based partly on Forest Inventory and Analysis (FIA) data. MODIS and MISR retrievals of forest fractional cover and mean height both showed compatibility with the USFS maps, with MODIS mean absolute errors (MAE) of 0.09 and 8.4 m respectively, compared with MISR MAE of 0.10 and 2.2 m, respectively. The respective MAE for aboveground woody biomass was ~10 Mg ha−1, the same as that from MISR, although the MODIS retrievals showed a much weaker correlation, noting that these statistics do not represent evaluation with respect to ground survey data. Good height retrieval accuracies with respect to averages from high resolution discrete return lidar data and matches between mean crown aspect ratio and mean crown radius maps and known vegetation type distributions both support the contention that the GO model results are not spurious when adjusted against MISR bidirectional reflectance factor data. These results highlight an alternative to empirical methods for the exploitation of moderate resolution remote sensing data in the mapping of woody plant canopies and assessment of woody biomass loss and recovery from disturbance in the southwestern United States and in parts of the world where similar environmental conditions prevail